Geometric Phase and Quantum Phase Transition : Two-Band Model 
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The connection between the geometric phase and quantum phase transition has been discussed 
extensively in the two-band model. By introducing the twist operator, the geometric phase can 
be defined by calculating its ground-state expectation value. In contrast to the previous numerical 
examinations, our discussion presents an exact calculation for the determination of the geometric 
phase. Through two representative examples, our calculation shows the intimate connection between 
the geometric phase and phase transition: different behaviors of the geometric phase can be identified 
in this paper, which are directly related to the energy gap above the ground state. 

PACS numbers: 75.10.Pq, 03.65. Vf, 05.30.Pr, 42.50.Vk 



I. INTRODUCTION 

Recently, quantum phase transition [l[ has received 
great attention due to its intimate correlation to the fun- 
damental principles of quantum mechanics .especially to 
the concept of quantum entanglement @, H, 0, [j| ( see 
Ref. @ for a comprehensive review ) . In general the quan- 
tum phase transition happens when the degeneracy of 
the ground states occurs, which cannot be characterized 
completely by the pattern of symmetry broken ( order pa- 
rameters of some kind ). Instead the universal quantum 
order or topological order is needed for the description of 
properties of the ground state in many-body systems Q • 
Recently the connection between the geometric phase of 
the ground state and the quantum criticality has been 
displayed in spin-chain systems by displaying the singu- 
larity of the geometric phase closed to the critical points 
[1, [9|. Moreover in Ref 0] the author showed that the 
scaling behavior of the geometric phase of the ground 
state near the critical points can also display the univer- 
sal class of the phase transitions. Many works have been 
devoted to this interesting issue E3, EH, E2, EH ( or see 
Ref. EH for a review ). 

Although great progresses has been made in the un- 
derstanding of quantum phase transition from the fun- 
damental principles of quantum mechanics, an impor- 
tant question is not yet resolved: whether there exists 
a universal way to characterize the different phases and 
their boundaries. More specifically, could the geometric 
phase of the ground state in many-body systems serve 
this purpose? This conjecture is natural since the quan- 
tum phase transition generally emerges from the degen- 
eracy of the ground state in many-body systems. Geo- 
metric phase as a measurement of the curvature of the 
Hilbert space, could mark the fantastic changes of the 
ground state when degeneracy happens. However the 



critical point is how to obtain the geometric phase of the 
ground state. The earlier method is to impose a local ro- 
tation about some special orientations, as has been done 
in [1, E3j- I n our own viewpoints, the connection built 
by this method is fragile; the geometric phase is trivial 
when the system is symmetrical about this rotation. 

Recently the differential information-geometry analy- 
sis of quantum fidelity in many-body systems displayed 
the intimate correlation between quantum phase tran- 
sition and the singularity of fidelity between the states 
across the transition point ■ In these papers the quan- 
tum geometric tensor, which is intimately related to the 
degeneracy of the ground state, was introduced for the 
determination of the fidelity. As shown in these papers, 
the imaginary part of this tensor actually described the 
curvature two-form whose holonomy is the Berry phase, 
and the degeneracy of the ground state would induce the 
singularity of fidelity E3- However, in our opinion, it 
seems not transparent to directly define the geometric 
phase in this coupling-parameter space since a cyclic evo- 
lution may be difficult to construct. Moreover, the ex- 
plicit expression for the ground state is necessary for the 
construction of this tensor, which in general is difficult. 
Furthermore it is also unknown in the case where the de- 
generacy of the ground state is broken. In Ref. EHj the 
Bargmann phase, a generalization of the Berry phase, has 
also been constructed for detecting the phase transition 
in many-body systems. However the connection between 
degeneracy of the ground state and the Bargmann phase 
was unclear in this case since there was a lack of simple 
interpretational underpinnings for the Barg mann phase 
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in terms of physical adiabatic processes 

With respect to the points stated above, it is urgent to 
find a popular way for construction of geometric phase 
in many-body systems. For this purpose, a nonlocal op- 
erator -the twist operator- is introduced to obtain the 
geometric phase of the ground state in this paper. Our 
calculation shows that the geometric phase, decided by 
the ground-state expectation value of the twist operator, 
can serve as the quantity to distinguish different phases 
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and boundaries between them, as shown below. The gen- 
eral form for the twist operator can be written as in the 
lattice systems, 



2ni ^ N 



(i) 



in which N is the number of lattice sites, x denotes the 
coordinates of lattice sites and n x is generally related 
to the physical quantities located at site x, such as the 
total spin or charge, or particle number at site x and 
so on. The twist operator was first introduced by Lieb, 
Schultz and Mattis for the proof of gapless excitation in 
one-dimensional spin-1/2 chains [15|. Then Resta pointed 
out that its ground state expectation was direct to the 
Berry-phase theory of polarization in strongly correlated 
electron systems 16|. Moreover Aligia found that the 



ground-state expectation value of the twist operator Eq. 
(dJ allows one to discriminate conducting from noncon- 
ducting phase in the extended quantum systems [TtJ • 
The vanishing of the ground-state expectation value, i.e. 
7] = 0, has been shown the ability to detect the bound- 
aries of different valence-bond-solid phases in spin chains 
[l8| . However these studies were implemented in some 
special examples and a general discussion was absent. 
Moreover since the previous calculations were numerical 
or approximate, the details of rj adjacent to the phase 
transition points are unclear. Hence it is of great interest 
to find the exact expression for rj, even for special cases. 
Our paper serves this purpose, and the exact results can 
be found for a special case. 

It is of great interest to note that the twist operator 
T) actually creates a wave-like excitation since it rotates 
all the particles with a relative angle between the neigh- 
boring lattices 2tt/N 15]. Under the large N limit the 
ground state has an adiabatic variation, and its ground- 
state expectation value is exactly a geometric factor, of 
which an argument is the geometric phase. Applying 77 to 
the unique ground state, one obtains a low-lying excited 
state. The important quantity is the overlap between the 
ground state and this excited state, i.e. 



(sMff>) 



(2) 



in which \g) denotes the many-body ground state. In 
general z is a complex number and its argument is just 
the geometric phase, determined by 



lg = Arg[z] = i / 
Jo 



d0< fl (0)|^|ff(0)), 



(3) 



in which \g((j>)) = exp(* ^2 X xn x )\g). Since 7 g in fact 
came from the continued deformation of the boundary 
condition of systems [13] and then slightly related to the 
symmetry of the Hamiltonian, this construction of the ge- 
ometric phase is more popular than the previous method. 
An important character is that 7 9 is related to the cor- 
relation functions for the ground state, and numerical 
evaluation could be implemented efficiently [lj| [l?], EB] • 



It is an immediate speculation that z and "f g may be 
singular near the critical points, where the degeneracy of 
the ground state happens and the macroscopic proper- 
ties of the system have fantastic changes. However, our 
findings are more subtle; the exact calculations show an 
unexpected ability for 7 S or z to distinguish the different 
phases in many-body systems; one case is that z tends to 
be zero and then j g is ill-defined when one approaches the 
phase transition points, in which the degeneracy of the 
ground state happens. The other is that 7 g has different 
values for different phases and displays the singularity at 
the transition points, where no degeneracy happens. The 
physical reason, as shown in the following discussions, is 
directly related to the energy gap above the ground state. 

The paper is organized as follows. In Sec. II, the ex- 
act expression of z and the geometric phase 7 ff are pre- 
sented in the two-band model for a special case, in which 
the ground state is the filled Fermi sea. In Sec. Ill, 
two representative examples are provided for the demon- 
stration of this connection. One is the D-dimensional 
free-fermion model, in which there are quantum phase 
transitions originated from the ground-state-energy de- 
generacy. The other example is the Su-Schrieffer-Heeger 
( SSH ) model, in which the energy gap is non-vanishing 
at the quantum phase transition points and a topological 
order, defined by the geometric phase 7 g , provides a clear 
description of the phase diagram for this model. 



II. TWO-BAND MODEL 

Consider the one-dimensional (ID) translational in- 
variant Hamiltonian with two bands separated by a finite 

gap mug, 



(4) 



x,x f 



in which c x ^ = (c_)_,c_)x defines a pair of fermion 
annihilation (creation) operators for each site x,x — 
1, 2, . . . , N and the form of c x is decided completely by 
the Hamiltonian. Jt? x ,x' is a 2 x 2 matrix and its elements 
can be determined by the hermiticity of Eq. ^ , 



C 



D 

A- 



(5) 



x,x' 



in which A± >xx > = A± , x and B XjX > = C*, . Although 
our discussion is restricted to ID system in this section, 
it should point out that this situation can be easily gen- 
eralized to higher dimension systems. 

Without the loss of generality, it is conventional for 
the translational invariant system to impose the periodic 
boundary condition. In spite of the simplicity, the Hamil- 
tonian Eq. ([J} has a wide range of applications, such as 
the Bogoliubov- de Gcnnes Hamiltonian in superconduc- 
tivity, graphite systems [l|| . Applying the Fourier trans- 
formation and considering the periodic boundary con- 
dition, Ca; = l/^/N^2 k e lkx Ck in which k — 2nn/N with 
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n = 1, 2, . . . , N. Then the Hamiltonian in the momentum 
space can be written generally as H = J^k c k J^(k)c k . If 
one introduces the four- vector i2^(fe)(/i = 0, x, y, z), then 
the Hamiltonian can be rewritten as, 



m k 



(6) 



in which ao is a 2 x 2 unit matrix and cr^ (i = x, y, z) is the 
Pauli operators. Obviously Eq. ([6|) can be diagonalized 
by finding the eigenvectors u± of J2i Ri{k)oi — R(fc) • er, 
in which the vector R(fc) is similar to the Bloch vec- 
tor for the density operator in the 2x2 Hilbert space 
[2l| . and furthermore should also satisfy the relation 
(R x (k), R y (k)) 7^ (0, 0) since in this case one has a trivial 
geometric phase [22}, 



1 



v± 



R x {k)-iR y {k) 

y/2R(k)(R(k)T Rz{k)) V ±R( - k ) ~ R *( k ) 



(7) 



in which R{k) = |R(fc)| and the corresponding eigenval- 
ues are E± — Rq (k) ± R z (k) . Obviously there is a finite 
gap between the two bands since E + > E_ . The ground 
state is defined as the filled Fermi sea \g) = Y\ k fe |0)fc, 
and |0)fc is the vacuum state of 



t - J, 



in which /3 

C±fc. 

Now it is time to determine the geometric phase Eq. 
©, given by the ground-state expectation value of the 
twist operator Eq. ([2]). In this model the twist operator 
can be expressed explicitly as 



(8) 



in which n x — c x c x = c^ + c + + c_c_. Seemingly c x c x 
could define the particle number at the site x, however 
the physical meanings for it may be different for different 
systems, dependent on one's interests; for spin systems 
it may denote the total spin at site x, and for electron 
systems it may also denote the total charge number at 
site x. The geometric phases in the both situations have 
been defined respectively as the spin Berry phase and the 
charge Berry phase, which have extensive applications 
in determining the phase diagram in strongly correlated 
electron systems [23j . 

It is a crucial step to determine z. First with the pe- 
riodic boundary condition, one can rewrite r\ in the mo- 
ment space 



k 



Now, introduce the new fermion operators 



(9) 



(10) 



in which v± is defined in Eq. (O and both of f3±^ satisfy 
the anti-commutative relation. Then the ground state is 



defined as the filled Fermi sea \g) = Yik P- fc|0)fe- Substi- 
tute Eq. CLO} into Eq. © 

r, = ]\exp[-^(filM0k + 0i ik dk0-,k+0l ik dk0+,k)], 

k 

(11) 

in which /3 k = (/3_,fc, P+,k) T and 



M = 



K' K 

~K* K'* 



(12) 



where 



cos Ok 
K" 
K 



R z (k) Ry(k) 

, 7fe = arctan ■ 



R(k) 



R x {k) 



-ism —dkik 



—j-{dk0k + isinOkdkjk) 



(13) 



The last two terms in Eq. (fTTj) precludes the further 
exact calculations. An important case is that if R y (k) = 
—R y {—k), one then can properly choose v± so that the 
new Fermi operator (3± yk can be converted into each other 
by exchanging k *-* —k. Then the terms in Eq. (jlip . 

P- k®k(3-,k,P\. kdk(3+ : k can cancel each other. An exact 
result in this special case can be obtained for z 



(ffWs>=ni i - 



i) 



\K\ 

C 2 _ 



-X- 



in which 



A± = ±i^/\K'\ 2 + \K\' 



ci = -\k\- 



{\ ± -K'f 



-1)] 
(14) 

(15) 



and from this formula, the geometric phase "f g can also 
be obtained exactly. 

The exact determining of z provides the ability to de- 
tect the distinguished behaviors of "f g or z near the phase 
transition points. Although the exact results can be ob- 
tained only for this special case, some different connec- 
tions between the geometric phase j g or z and quantum 
phase transitions are disclosed, as shown in the following 
calculations. 



III. EXEMPLIFICATIONS 

In this section two representative models are presented 
to display the distinguished characteristics of "f g or z 
for the determination of the phase diagram in many- 
body systems. One is the D-dimensional frcc-fermion 
model, in which the quantum phase transition is orig- 
inated from the degeneracy of the ground-state energy 
[H]. The other is the Su-Schrieffer-Heeger (SSH) model, 
in which the quantum phase transition happens with the 
non- vanishing energy gap above the ground state [l9l |. 
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It is obvious that these examples include two important 
cases of quantum phase transition; one is the degeneracy 
of the ground-state energy, whereas not for the other. 
Moreover the both models are exactly solvable. 



D-dimensional free-fermion model 



The Hamiltonian is read as 



H 



<iJ> 



7(4 C ] 



h.c.)] 



(16) 



in which denotes the nearest-neighbor lattice sites 
and Cj is the fermion operator. This Hamiltonian, first 
introduced in Ref.[24li depicts the hopping and pairing 
between nearest-neighbor lattice sites, in which A is the 
chemical potential and 7 is the pairing potential. Eq. (|16[) 
could be considered as a D-dimensional generalization of 
one-dimensional spin-1/2 XY model. However for the 
D > 1 case, this model shows some novel phase charac- 
teristics |24j |. 

Eq. p6p can be resolved exactly by transforming into 
moment space with periodic boundary condition. Wiht 
the help of the Bogoliubov transformation [24| , one has 



const. 



(17) 



in which Ak = y/t 



f k = E„=i cosfc « - A and 
Ak = 7X)a=i sm ^a- The phase diagram can be de- 
termined based on the gapless excitation = [24j . 
For D = 1, which corresponds to the one dimensional 
spin-1/2 XY model, the energy gap above the ground 
state is non-vanishing except at A c = 1 for 7 / 0, where 
a second-order quantum phase transition occurs. For 
7 = the energy of the ground state is degenerate in 
the region |A| < 1 and the transition occurs at A = ±1. 
When D = 2, the phases diagram should be identified 
with respect to two different situations; for 7 = 0, the 
degeneracy of the ground state occurs when A € [0,2], 
whereas the gap above the ground state is non- vanishing 
for A > 2. However for 7 / three different phases can 
be identified as A = 0, A e (0, 2] and A > 2. The first 
two phases correspond to case that the energy gap for 
the ground state vanishes, whereas not for A > 2. One 
should note that 7 = means a well-defined Fermi sur- 
face with k x = k y ±TT, whose symmetry is lowered by the 
presence of 7 terms. For D = 3 two phases can be iden- 
tified as A G [0, 3] with the vanishing energy gap above 
the ground state and A > 3 with a non-vanishing energy 
gap above ground state. In a word the critical points can 
be identified as A c = D(D = 1,2,3) for any anisotropy 
of 7, and A = for D = 2 with 7 > 0. 

Defining the fermion-pair operator Ci = (c^c^f and 
transforming the system into moment space, one then 
obtain 



H 



A - X) Q C0S k ° 



-A + J^a C0S k o 



c k (18) 
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FIG. 1: ( Color online ) The ground-state expectation value 
of twist operator 2 for the D-dimensional free-fermion model 
vs the parameter A. We have chosen 7=1 for this plot. The 
plotting in the inset of (c) shows the details of A > 3. 
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It is worth noting R(k) = (0, 7^ a sinfc Q , A — J2 a cos k a ) 
and R(k) — A^. R y (k) is obviously satisfied the require- 
ment i? y (k) = -Ry(-k). Substituted into Eq. (fT4| . one 
obtain 



D 



lg = Argz M Arg J [ J]_ cos — | K t 



a — 1 kr, ,k,. 



N 



(19) 



in which /i = 1 ... D and 



~ J2 a cosfc a )[cosfc M (A - J2 a cos k a ) - sinfc M ^ sinfc a ] 
2[ 7 2(^ Q sin fc Q )2 + (A - E Q cos fc Q ) 2 ]3/ 2 

I 



(20) 
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FIG. 2: ( Color online )(a) z for A = 1 ( plotted by a logarithm 
) vs the particle number N for D = 1; (b) the asymptotic 
behavior of z closed to transition point A = 1 for D = 1; (c) 
the geometric phase 7 9 for D — 1. We have chosen 7 = 1 for 
all plots. 



The schematic drawing of z for D = 1,2,3 have been 
presented in Figs. [TJ in which we have chosen 7=1 for 
specification. Some different characters can be found in 
the figures [25j . 

D = 1. It is the well-known one-dimensional XY 
model for this case, in which there is a quantum phase 
transition at A = 1 because of the degeneracy of the 
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FIG. 3: z for D = 2 vs particle number TV at phase transition 
point A = 2. 



ground-state energy, z are plotted with A respectively 
in Fig. [T] (a). It is obvious that z has an dropping and 
then an abrupt increment when one approaches the phase 
transition point A c = 1. Moreover our calculation shows 
that z tends to be zero with A — > 1 as shown in Fig. [2] 
(b) , and at exact A c = 1 z tends to be 1 with the increase 
of the lattice site number as shown in Fig. [5] (a) [26[ . j g 
have also been plotted in Fig. [2] (c), in which a rapid 
oscillation happens closed to A c = 1. These phenomena 
mean that j g is ill-defined closed to A = 1 and has an 
abrupt change at the phase transition point. Since the 
energy gap vanished only at A c = 1 , the singularity of j g 
and z would be directly related to the degeneracy of the 
ground state energy. It also hints that one could mark 
the transition point by detecting the point where 2 = 0. 

D = 2. With the increment of dimensionality, the 
situation becomes more complex. We have plotted z in 
Fig. [1] (b) . It is obvious that two different regions can 
be identified as A € [0, 2] in which z is disordered, and 
A > 2 in which z is an increasing function of A. With re- 
spect that the disappearance of the energy gap above the 
ground state happens when A S [0, 2], z presents a clear 
identification of the phase diagram. It is a reasonable 
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speculation from Fig. Q] (b) that z may tend to be zero 
under the thermodynamic limit when A E [0, 2]. Then j g 
is under thermodynamic limit 



la 



\ undetermined, A < 2 
|0, A>2. 



(21) 



Our calculation also shows that z tends to be zero with 
the increment of N at the exact transition point A c = 2, 
as shown in Fig. [3] Similar to the case of D = 1, it may 
be desirable to find the point z = as a way of detecting 
the phase transition. 

D = 3. This case is very similar to that of D = 2, 
except the phase transition happens at A c = 3. z has 
been shown in Fig. [1] (c). However in this case z seems 
unlikely to detect the phase transition since the data in 
the figure has a smoothing changes at the phase transi- 
tion point for large N, as shown in the inset of Fig. [T] 
(c). Only for N = ll 3 , there is an abrupt changs of z 
near to the phase transition point. One should note that 
z tends to be zero when A G [0,3], in which the energy 
gap above the ground state disappears. 

From the discussions above, one can note the great im- 
pact of the degeneracy of the ground-state energy on the 
geometric phase j g or z; The degeneracy of the ground- 
state energy leads to z = or the ill-defined j g . However, 
this conclusion would not be made until the next exam- 
ple is studied in which the energy gap above the ground 
state does not disappear. It is interesting to give a fur- 
ther discussion of the geometric phase in this nontrivial 
case. 

Unfortunately Eq. (|20[) seems unsuccessful in char- 
acterizing the transitions for 7 = (the tight-binding 
model) since in this case is completely undetermined 
when A — ^ Q cos k a = 0. With respect to Eqs. (jTHJ) and 
flS]), it means that [H, rj\ = since E(ij> c i c j; J2i c i c i] = 
in this special case. Hence our discussion excludes this 
special situation since one has trivial results. One should 
note that the phase of A = with 7 > in D = 2 also 
cannot be identified by z since the transition comes from 
the deformation of the Fermi surface instead of the de- 
generacy of the ground state. 



B. Su-Schrieffer-Heeger (SSH) model 

Another example is the Su-Schrieffer-Heeger (SSH) 
model, which is also exactly solvable. The ID tight- 
binding Hamiltonian for the SSH model for a chain of 
polyacetylene is given by [27[ 



N 



(22) 



2=1 



in which <fii represents the dimerization at the Ith site 
and an alternating sign of the hopping elements reflects 
dimerization between the carbon atoms in the molecule. 



1- 


-o-N=11 <&A 




N=101 W 




-□-N=1001 \ 




0- 


















-1- 











1.0 



0.5 



-o-N=11 
-a-N=101 
-□-N=1001 





FIG. 4: ( Color online 
parameter <f>. 







7 9 and zj for the SSH model vs. the 



Without the loss of generality, it is convenient to neglect 
the kinetic energy in the system and take cf>i = <j), t = 1 
[20I ] . There is a critical point at 4> — , which divides 
the ground states into two different types. One should 
point out that SSH has a gaped excitation for any ip G 
[— 1 , 1] and the quantum phase transition comes from the 
excitation of the boundary states [20]. 

One can find from the following calculation that j g 
can discriminate the two phases and the boundary be- 
tween them. Defining c x — (c;,c;+i) T , the Hamiltonian 
becomes 



H 



N 

£4 

x=l 

t 



-(1 



-(1 



Cz+i + h.c. 



(23) 



Imposing the periodic boundary condition and Fourier 



7 



transformation, we then have 



Jf(k) 



_(i + 0_(i 

-(l + 0)-(l-<£)e ife 



-, — ik 



(24) 

It is obvious that R(fc) = (— (l + (f>) — (1 — <t>) cosfc, — (1 — 
</>)sinfc,0) and R y (k) = —R y (—k) is satisfied. Then 

-pp. 2V2tt zlm[iC] . Sygrr 

cos \K \ = sin \K ), (25) 



in whirl, - 1 (l-0) 2 + (l-0 2 ) cos fc 

We plot 7 g against with different site numbers in 
Fig. [4j It is obvious that -f g is — n for <f) £ [— 1, 0) and 
zero for </> G (0, 1] from Fig. [5] Moreover 7 S tends to be 
7r at exact phase transition point = 0. Furthermore 
our calculation shows that \z\ is not zero, which means 
that one cannot detect the phase transition by finding the 
point z = 0. Since the energy gap for the ground state is 
nonvanishing in this model, the geometric phase j g can 
be well-defined for any <j>. With respect to the discussion 
for the /^-dimensional free-fermion system, it is evident 
that 7 S or z are directly related to the degeneracy of the 
ground state. 



IV. DISCUSSIONS AND CONCLUSIONS 

Given the two examples, some comments should be 
presented in this section. In this paper the twist oper- 
ator Eq. |T]) has been introduced and its ground-state 
expectation z has been calculated to define the geomet- 
ric phase Eq. © for the two-band model. Although the 
absent of general results, the exact expression of z can 
be obtained in a special case, which provides the abil- 
ity to detect the details of the geometric phase adjacent 
to the phase transition points. With respect to the dis- 
cussions for two representative examples-D-dimensional 
free-fermions model and the Su-Schrieffer-Heeger model, 
some distinguished properties of 7 9 or z have been found 
in our calculations. 

First when the degeneracy of the ground state hap- 
pens, z tends to be zero and the geometric phase j g is 
ill-defined in this shown in Figs. [TJ while j g can 

be well-defined when the energy of the ground state is 
nondegenerate, as shown in Fig. 2] This phenomenon 
clearly displays the intimation connection between ge- 
ometric phase, defined by the ground-state expectation 



value of the twist operator, and the degeneracy of the 
ground state in many-body systems. Consequently one 
can find the nodal structure of the geometric phase (the 
situation that the geometric phase is ill-defined because 
of z = 0) to detect the phase transition originated from 
the degeneracy of the ground state. The nodal structure 
of geometric phase is introduced by Fillip and Sjoqvist 
for the description of experimental measure of the geo- 
metric phase based on the interference, which character- 
izes the condition for the disappearance of the fringes and 
then the geometric phase is ill-defined. Second geometric 
phase can also present the phase diagram even if there is a 
energy gap above the ground state, as shown in Fig. [4J in 
which two different phases defined by 7 g = — 7r and can 
be identified and the phase transition point is marked by 
the discontinued variation of geometric phase. In a word 
the geometric phase 7 ff displays the ability to mark the 
phase diagram in this discussion, whether the phase is de- 
termined by the degeneracy of the ground state or not. 
Hence the geometric phase can provide a more popular 
depiction for the phase transition. 

However there still exist some problems. First is that 
the geometric phase seems to fail to characterize the 
tight-bond model (7 = in Eq. (HU)). It is the rea- 
son that [if, 77] = 0, and the twist operator has a trivial 
effect on the ground state. Second the geometric phase 
fails to detect some phase transitions not originated from 
the degeneracy of the ground state, such as the transi- 
tion from the deformation of the Fermi surface. Thirdly 
7 g or z seems unable to detect the broken of symmetry 
which happens in the ID spin-1/2 XY model, as shown 
in Fig. [2(a). Although there exists some defects, the ge- 
ometric phase defined by the twist operator provides one 
another way of detecting the phase diagram for many- 
body systems. Moreover the flexibility of choosing n x 
in the definition of the twist operator implies that one 
could properly choose different physical quantities n x for 
the description of different properties of the system. 

Note added. Recently we become aware of a paper 
which also focuses on the connection between the geomet- 
ric phase and the quantum phase transition by numerical 
evaluation [29]. In this paper, three different phases in 
gapped spin chains can be defined by the geometric phase 
7 = 0, 7r, and undefined respectively, which is similar to 
our conclusions. 
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